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The Minimum Variability Time Scale and its Relation to 
Pulse Profiles of Fermi GRBs 
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ABSTRACT 

We present a direct link between the minimum variability time scales extracted 
through a wavelet decomposition and the rise times of the shortest pulses extracted 
via fits of 34 Fermi GBM GRB light curves comprised of 379 pulses. Pulses used in 
this study were fitted with log-normal functions whereas the wavelet technique used 
employs a multiresolution analysis that does not rely on identifying distinct pulses. By 
applying a corrective filter to published data fitted with pulses we demonstrate agree- 
ment between these two independent techniques and offer a method for distinguishing 
signal from noise. 
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1 INTRODUCTION 



One approach for probing light curves whic h has 
received attention jNemiroflJ l200fj|; [Norris et alj 120051; 



IHakkila fc Nemirofdl2009l : lNemirofdl2012T ) is to express them 
as a series of displaced pulses, each with a parametric form. 
There is an appeal to this approach because fitting routines 
are well-understood and interpretations of rise time, decay 
time, full width at half max, etc, are possible. On the other 
hand, one must make certain assumptions when using the 
pulse-fitting procedure such as the choice of the functional 
form to use for an individual pulse and the number of pa- 
rameters to be included in the fitting function. Moreover, 
light with high variability at low power may show variations 
which are not statistically significant. 

A complementary approach using a wavelet-based anal- 
ysis of a set of both long and shor t GRB light curves is dis- 
cussed bv lMacLachlan et alj ((2012) in which a time scale, Tp, 
is identified that marks the transition from white noise to a 
power law in the power density spectrum (a f~ a behavior), 
ft is argued that over time scales smaller than rg the light 
curves appear stochastic and signal power is distributed uni- 
formly. At time scales larger than rg, identifiable structures 
(such as pulses) with signal power are no longer distributed 
uniformly over the periods of light variation. For this reason 
Tp is referred to as the minimu m variability time scale. 

The analysis presented in (|MacLachlan et all I2012T I is 
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a non-parametric approach to probing light curves for time 
scales. It makes no assumptions about the nature of the 
structures in a given light curve that give rise to the f~ a 
character. The technique, however, offers no firm connec- 
tion between rg and the constituent structures although it 
seems reasonable to associate rg with the scale of the small- 
est emitting structures present. 

Results from an application of a log-normal pulse- 
fit ting procedure to GRB light curves have been published 
bv lBhat et all (|2012f l. In this paper we make a meta-an alysis 
of the timing results presen ted bvlBhat et alj (|2012|) com- 
pared with the techniques of IMacLachlan et alj ( 20121 ) for a 
set of 34 GRBs used in both studies. 



2 ANALYSIS 

We begin by considering the relati on between rg an d the 
pulse parameters given in Table 3 o f iBhat et all (|2012l ). The 
parameters with temporal units in Table 3 are: time-since- 
trigger, rise time, decay time, and FWHM. In all, 34 GRBs 
comprising 379 pulses are considered here. We note that 
ri se time, decay tim e, and FWHM as presented in Table 3 
of lBhat et all (|2012l ) are tightly correlated and for the argu- 
ment that follows are interchangeable without affecting the 
conclusions . However, we use rise time to make our argument 
because, as IBhat et alj l|2012T ) noted, rise times are observed 
to be short e r tha n decay time and FWHM (see Table 3 in 
IBhat et alj (|2012r >). We considered only those light curves 
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Figure 1. Minimum variability time scale versus rise time in 
th e Observer frame . The rise times are taken from Table 3 
of iBhat et al. | j2012h . The line represents the locus where r« = 
rise time. Wc identify the area below the line with white noise. 
The data are expected to press up against the line from above 
but not to cross it. 



from Nal detec tors and summed o ver the energy acceptance 
as m Table 3 of lBhat et all (j2012T ) and in lMacLachlan et al.l 
(2012h . 

In Fig. [1] we plot the rise time for all 379 pulses (34 
GRBs) along the vertical axis and Tp along the horizontal. 
Note that for each GRB for which one rg is computed, there 
is a possibility for multiple pulses and therefore multiple rise 
times, hence the vertical columns of rise times for a single 
value of Tp. For a given column of pulse times the short- 
est pulse rise times are at the bottom and one finds larger 
rise times by moving up the column. An equality line is also 
shown which is the locus where 773 equals rise time. Arguing 
as we do that Tp represents the minimum variability time 
scale the space in the r^-rise time plane below the equality 
line should be interpreted as a structureless white noise re- 
gion. If some method were capable of discerning light curve 
structure in the region we define as white noise, then our as- 
sertion of having found a minimum variability time scale will 
have been disproven. Indeed, in Fig. [T] there are 27 pulses 
with rise times below the equality line. The uncertainties 
accompanying these 27 rise times are small, making their 

intrusion into the white noise region significant. 

However, a closer inspection of Table 3 of IBhat et al.l 

(|2012l ) reveals that there are 20 light pulses in Fig.[T]with rise 
times that are smaller than the smallest bin widths, in some 
cases smaller by factors of ten or a hundred. Moreover, of 
those 20 pulses there are 16 pulses in Fig. [T] with full widths 
at half max that are smaller than the smallest bin widths 
and indeed those 16 all fall below the equality line. While it 
seems that inclusion of these pulses in Table 3 is important 
for the sake of completeness, w e question the phys i cal rea lity 
of these pulses. Note that in iMacLachlan et all |2012T ) all 
light curves are binned at 200 microseconds. Fig. [2] shows 
the effect of removing the 20 non-physical pulses. Note that 
in Fig. [2] the white noise region has been vacated by all 
but seven points and none of the pulse rise times above the 
equality line have been disturbed by the bin width cut. For 
the seven points that remain beneath the equality line, we 
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Figure 2. Minimum variability time scale versus rise time in the 
Observer frame. We have removed all pulses (20) with rise times 
smaller than the light curve bin width. 



Rise Time vs x. 



100 



























' White Noise 



0.001 0.01 0.1 1 10 

Xp Observer Frame [sec] 

Figure 3. Minimum variability time scale versus rise time in the 
Observer frame as in Fig. [2] We have folded a single bin width 
into the rise time uncertainties. 

show in Table[T]that six are within one sigma of the equality 
line. 

We make one other point regarding the pulse rise times 
in Fig. [1] and Fig. [2] in particular regarding the size of 
th e uncertaint i es. O f the 379 pulse rise times reported 
by IBhat et al.l (|2012h and used for this meta-analysis, 301 
have uncertainties smaller than the binning of the light 
curve, in some cases hundreds or thousands of times smaller. 
We argue that a conservative estimate of the uncertainties 
for the pulse rise times should be no smaller than a bin 
width. Thus, we add in quadra ture a bin width (as reported 
in Table 3 of lBhat et al] (|2012T ) ) to the rise time uncertain- 
ties (also reported in Table 3 of lBhat et al.l l|2012l )) and plot 
the result in Fig. [3] 

In Fig. [4] we plot only the smallest rise times for each 
GRB against 773. We argue that by rejecting pulse rise times 
smaller than light curve bin widths and by folding rise 
time uncertainti e s with a bi n width we get strong ev idence 
that IBhat et al] (|2012l ) and IMacLachlan et al] (|2012T ) have 
tracked the same physical observables over approximately 
three orders of magnitude using independent methods. 
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Table 1. Pulses with rise times smaller t han Tp but larger than bin widths. Pulse Number (#), rise (time), tSrise (time), and bin width are 
taken from Table 3 of lBhat et al. The column labeled Arise (time) is obtained by combining <5rise and bin width in quadrature. 

The columns diff and % diff refer to the differences between Tp and rise (time). 



GRB Pulse # Tp [s] St~ [s] 8tJ [s] rise [s] <5risc [s] Arise [s] bin width [s] diff [s] % diff 
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Figure 4. Minimum variability time scale versus rise time in the 
Observer frame as in Fig. [3] but with only smallest rise times 
included. Note that the equality line between Tp and rise time 
marks a boundary between scaling processes and white noise and 
gives substance to the interpretation of the minimum variability 
time scale. 



3 RESULTS AND DISCUSSION 

For a large sample of s hort and long Fermi GBM bursts, 
iMacLachlan et al] l|2012l ) used a technique based on wavelets 
to determine the minimum variability time scale, Tp. The 
authors associate this time scale with a transition from red- 
noise processes to parts of the power spectrum dominated by 
white noise or random noise components. Accordingly, the 
authors note that this time scale is the shortest resolvable 
variability time for physical processes intrinsic to the GRB. 
In addition, histograms of the values of Tp for long and short 
GRBs were shown to exhibit a clear temporal offset in the 
mean Tp values for long and short GRBs. 

In a separate analysis, using a p articular functional form 
for pulse shapes, iBhat et al] (|20T2t ) have extracted an ex- 
tensive set of key pulse parameters such as rise times, decay 
times, widths (FWHM), and times since trigger for a host 
of bright GRBs detected by Fermi/GBM. Using the FWHM 
values, these authors also reported a significant temporal 
offset between the mean values for long and short GRBs. 

Although neither group offers a concrete explanation 
for the temporal difference between the distributions of long 
and short GRBs, it is noteworthy that they arrive at a result 
which is quantitatively in good agreement with one another, 
especially having used independent approaches. Both sets of 



analyses also suggest scaling trends between characteristic 
timescales. In the case of minimum variability timescales the 
trend is between Tp and the duration of the burst, typically 
denoted by Tgo. For the pulse-shape analysis, the trend is 
more readily evident and is demonstrated through a number 
of positive correlations involving key parameters such as rise 
times, decay times and FWHM times. 

It is relatively straightforward to interpret the scaling 
trends in terms of the internal shock model in which the 
basic units of emission are assumed to be pulses that are 
produced via the collision of relativistic shells emitted by the 
central engine. In the case of the pulse-fittin g method this is 
essent ially the default assumption. Indeed, iQuilligan et al] 
( 2002 ) in their study of the brightest BATSE bursts with 
T90 > 2 sec were the first to demonstrate this explicitly by 
identifying and fitting distinct pulses and showing a strong 
positive correlation between the number o f pulses and the 
durat io n of the burst. Mor e recent studies, Bhat k, Guiried 
(|201l] ); lHakkila fc Cumbeel (|200Sl ); lHakkiIa fc Preecd (|201ll ) 



have provided further evidence for the pulse paradigm view 
of the prompt emission in GRBs. 

The wavelet analysis does not, however, rely on iden- 
tifying distinct pulses but instead uses the multiresolution 
capacity of the wavelet technique to resolve the smallest 
tempor al scale present in the p rompt emission. Nonethe- 
less, as IMacLachlan et al] (|2012t ) have demonstrated, if the 
smallest temporal scale is due to pulse emissions, then we 
can still get a measure of the upper bound on the number of 
pulses in a given burst through the ratio T90 jrp . In the sim- 
ple model in which a pulse is produced every time two shells 
collide, the ratio Tqo/tp, should show a correlation with the 
du ration of the burst . Indee d, this correlation was reported 
by IMacLachlan et al] (|2012h . 

The similar trends of scaling demonstrated by these two 
methods, not only suggest the robustness of both methods, 
but also point, perhaps more importantly, to an underly- 
ing interconnection between key parameters extracted by 
the two techniques. In other words, the minimum variabil- 
ity time scale extracted by the wavelet technique is directly 
related to key pulse time parameters such as rise times (as 
depicted in Fig. under suitably controlled pulse-fitting 
methods. 



4 CONCLUSIONS 



Throu gh a met a-analysis of results prese nted bv lBhat et al] 
J2012I ) and by IMacLachlan et al] (|2012l) . we have studied 
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the relationship between key parameters that describe the 
temporal properties of a sample of prompt-emission light 
curves for long and short-duration GRBs detected by the 
Fermi/GBM mission. We compare the minimum variability 
timescale extracted through a technique based on wavelets, 
with the pulse-time parameters extracted through a pulse- 
fitting procedure. Our main results are summarized as fol- 
lows: 

a) Both methods indicate a temporal offset between 
short and long-duration bursts. The quantitative agreement 
between the two methods is quite good. 

b) Both methods point to scaling trends between char- 
acteristic timescales. In the case of the pulse-fitting method 
the scaling appears to involve parameters such as rise times, 
FWHM, and pulse intervals. For the wavelet technique, the 
scaling involves a correlation between the minimum variabil- 
ity time scale and the duration of the bursts. 

c) By demonstrating a strong positive correlation be- 
tween t/3 and the rise time of the shortest fitted pulses, we 
provide for the first time, a direct link between the short- 
est resolvable temporal structure in a GRB light curve with 
that of a key pulse profile parameter. 

d) By combining the two techniques, we have shown 
that one can arrive at a much tighter demarcation of the 
boundary between the power spectrum domains that sepa- 
rate red noise and white noise processes. 



5 ACKNOWLEDGEMENTS 

The NASA grant NNX11AE36G provided partial support 
for this work and is gratefully acknowledged. The authors, 
in particular GAM and KSD, acknowledge very useful dis- 
cussions with Jon Hakkila and Narayan Bhat. 



REFERENCES 

Bhat, P. N., & Guiriec, S. 2011, BASI, 39, 471 
Bhat, P. N. et al. 2012, ApJ., 744, 141 (online- version) 
Hakkila, J. & Cumbee, R. 2008. larXiv:0901.3171V l [astro- 
ph.HE]. 

Hakkila, J. & Nemiroff, R. J. 2009, ApJ., 705, 372 
Hakkila, J. & Preese, R. 2011, ApJ., 740, 104 
MacLachlan, G. A., et al. 2012. larXiv: 1201. 4431 V 2 [astro- 
ph.HE] 

Nemiroff, R. J. et al. 2000, ApJ., 544, 805 
Nemiroff, R. J. 2012, MNRAS, 419, 1650 
Norris, J. P. et al. 2005, ApJ., 627, 324 
Quilligan, F., et al. 2002, A.& A., 385, 377 



